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Abstract 

We study the linearised stability of the nakedly singular negative mass Schwarzschild 
solution against gravitational perturbations. There is a one parameter family of pos- 
sible boundary conditions at the singularity. We give a precise criterion for stability 
depending on the boundary condition. We show that one particular boundary condi- 
tion is physically preferred and show that the spacetime is stable with this boundary 
condition. 
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1 Introduction: Phantoms and Runaways 



The positive mass theorems are generally taken as a triumphs of modern relativity theory. 
They establish that, under plausible physical assumptions, asymptotically flat solutions of 
the Einstein equations with physically acceptable matter sources cannot have negative total 
mass. There is no shortage of solutions with negative mass, for example the Schwarzschild 
solution with negative mass parameter, but they are typical nakedly singular and hence 
assumed to be physically unacceptable. They are not expected to arise from regular initial 
conditions. However, this begs the question of why exactly it is that we are reluctant 
to consider solutions with negative mass. The question acquires topicality from various 
suggestions by cosmologists that the observed acceleration of the scale factor of the universe 
may be due to 'phantom matter', that is some matter, typically a scalar field, with negative 
kinetic energy [1, 2]. 

The idea of anti-gravity is, of course, a staple of science fiction writers. According to 
Mach [3] it appears to have been Foppl [4] who first, by analogy with electrostatics, explored 
the idea that gravitational masses could be both positive and negative. The discovery of 
dark energy may be said to establish anti-gravity as a serious subject for scientific discussion 
and the Randall-Sundrum scenario I [5], with its negative tension branes, has only reinforced 
the trend of considering matter with exotic energy momentum tensors. 

One problem with negative masses was pointed out long ago in a beautiful analysis 
of Bondi [6]. He drew attention to some special features of negative masses in general 
relativity. 

Firstly, as a consequence of the weak equivalence principle, a particle of negative mass 
falling in a gravitational field should fall at the same rate and in the same direction as a 
particle of positive mass. At the level of the geodesic equations this is because the mass 
cancels out from the equations of motion. Thus, for example, a cloud of negative mass 
particles, let's call them ghosts or phantoms, should accrete onto ordinary black hole of 
positive mass just like ordinary particles. 

On the other hand, a body with negative mass should repel particles with either negative 
or positive mass. Again, at the level of geodesic equations this is easily verified for the 
negative mass Schwarzschild solution. Combining this fact with our first observation, we 
see that the interesting possibility of a runaway solution exists in which a positive mass 
particle is chased in some direction by a negative mass particle, the combined system going 
into a state of constant acceleration. The positive mass particle attracts the negative mass 
particle to itself, but at the same time the negative mass particle pushes the positive mass 
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particle away. See Figure 1. 
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Figure 1. A negative mass particle chasing a positive mass particle. 

Bondi showed that the runaway phenomenon actually arises by exhibiting exact runaway 
solutions of the Einstein equations. His paper was written before the development of black 
hole theory, but may readily be adapted to include black holes [7]. Bondi considered the 
axially symmetric static vacuum metrics first studied by Weyl 



The vacuum equations are satisfied if the Newtonian potential U = U(p, z) is harmonic with 
respect to the standard Laplacian on three-dimensional Euclidean space E 3 with cylindrical 
polar coordinates (p,(p,z). The function k(p,z) may then be obtained by quadratures. A 
physically acceptable solution must have k = on any portion of the axis of symmetry to 
avoid conical singularities. This places restrictions on the possible solutions. The general 
formalism is reviewed in [8]. The Schwarzschild solution with positive mass M is obtained 
by taking for U the Newtonian potential of a uniform rod of length L = 2M. See Figure 
2. The portion of the line p = occupied by the rod is then a regular event horizon. The 
standard Schwarzschild coordinates are given by a system of confocal prolate ellipses in the 
p — z plane, with the rod a degenerate member of the family. 

A rod of mass per unit length | and of infinite length, occupying the positive z axis for 
example, corresponds to a Rindler horizon. Superposing a positive mass rod along the neg- 
ative z axis gives a solution with a conical defect somewhere on the axis of symmetry. This 
solution is the well known C-metric. Again, see Figure 2. These same conical defects arose 
in the solutions considered by Bondi which had as sources a continuous matter distribution 
and which were meant to model stars. Bondi's observation [6], which also applies to the 
black hole case [7], was that by superposing a further source in between the positive mass 
Schwarzschild rod and the Rindler rod the conical defects may be eliminated as long as the 
new source has negative mass. In this way a uniformly accelerating runaway solution may 
be constructed. The simplest negative mass solution for U (p, z) to take is that of a negative 



ds 2 = _ e 2U dt 2 + e -2U e 2k (dp 2 + + p 2 # 2 



(1) 
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mass point particle which results in a singular solution called the Curzon solution, but there 
are many other possibilities. In particular, one could use the negative mass Schwarzschild 
solution. This would amount to taking the Newtonian potential of a uniform rod of mass 
per unit length —5. This configuration is shown in Figure 2. 
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Figure 2. Rods at p = corresponding to: (a) A positive mass Schwarzschild black hole, 
(b) a Rindler horizon, (c) an accelerating black hole with a conical defect [the C-metric], (d) 
a negative mass black hole chasing a positive mass black hole in a nonsingular spacetime. 

It seems clear from the discussion above that under suitable circumstances the negative 
mass Schwarzschild solution should be unstable. As derived above the instability is dynam- 
ical and rather nonlinear. However, it is an interesting question whether this instability, 
or indeed other possibly unrelated instabilities, appear in a linearised analysis of classical 
perturbations around the Schwarzschild phantom. As far as we are aware, a discussion of 
the linearised stability of the negative mass Schwarzschild solution against gravitational 
perturbations has not been given before, and the purpose of the present article is to provide 
one. 

The crucial subtlety in our analysis will concern boundary conditions at the singular- 
ity. The linearised Einstein equations reduce to a set of Schrbdinger-like equations with a 
time-independent 'Hamiltonian'. The choice of boundary conditions is constrained by the 
requirement of a self-adjoint Hamiltonian, and hence a unitary time evolution. For scalar 
metric perturbations there is not a unique boundary condition 1 . We shall find a precise 
range of boundary conditions that guarantee linearised stability of naked singularity. The 



1 The problem of boundary conditions in negative mass Schwarzschild spacetime has been studied for test 
scalar fields in [9, 10, 11]. 
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existence of boundary conditions giving linearised stablility of the negative mass singularity 
is somewhat counterintuitive. Our results have a very similar flavour to recent results on 
the stability of Anti-de Sitter space, which is also not globally hyperbolic and for which an 
ambiguity of boundary conditions exists at infinity [12]. 

In the Randall-Sundrum scenario I [5], the negative tension brane is stabilised against 
increasing its area by a Z2 quotient of spacetime that fixes the brane. A similar phenomenon 
stabilises orientifolds in string theory and their supergravity realisation in terms of the 
Atiyah-Hitchin metric. The linearised stability of negative mass black holes that we have 
found here does not depend on any such discrete quotient. 

The organisation of this paper is as follows. We review the Schrodinger equations de- 
scribing perturbations of the spacetime. There is seen to be a one dimensional family of 
possible boundary conditions at the singularity We derive a critical boundary condition 
separating stable from unstable spacetimes. There is one particular boundary condition 
which results in perturbations with a more physical behaviour than the others. The space- 
time with this boundary condition is stable. We end with a discussion of our results. 



2 Stability analysis 

2.1 Negative mass Schwarzschild spacetime 

We will consider linearised perturbations about the four dimensional Schwarzschild space- 
time. We write the metric as 

ds 2 = -f{r)dt 2 + JJ^dr 2 + r 2 do\ 2) , (2) 

where da 2 2 ^ is the metric of a unit sphere and 

/(r) = l + £. (3) 

We are interested in the negative mass case, so fi = —2M > 0. The negative mass 
Schwarzschild spacetime is well known to have a naked singularity at r = 0. This is il- 
lustrated in Figure 3, which shows the Penrose diagram for the spacetime. 
As usual, it will be convenient to introduce the Regge- Wheeler coordinate 

dv I t \ 

As r — > 0, f(r) diverges. Therefore the range of the Regge- Wheeler coordinate is < < 
00. The timelike singularity is located at r* = and infinity corresponds to r* — > 00. 
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Figure 3. Penrose diagram for the negative mass Schwarzschild: There is a timelike 
singularity at r = 0. The boundary at r — > oo consists of two null lines. Anywhere in this 
region the Killing field d/dt is timelike. 

2.2 Gravitational perturbations 

Gravitational perturbations, g a t — > g a b + habi of the four dimensional background metric (2) 
may be grouped into two types: those of axial (vector) and polar (scalar) type with respect 
to their parity transformation on the sphere. Different types of perturbation do not mix at 
the linearised level. 

In order to solve the linearised Einstein equations, there is a standard technique for sep- 
arating the angular coordinates and constructing gauge invariant scalar variables which we 
will denote universally by $ [13, 14]. The quantities are linear combinations of perturbed 
metric components and their derivatives: h,d r h,d rr h. For each type of perturbation, one 
obtains an equation of motion for $ of the form 

_yV (5) 

The potentials appearing in such equations were first derived in a fixed gauged by Regge- 
Wheeler [15] and by Zerilli [16]. A unified treatment is given by Chandrasekhar [17, 18]. 
The gauge invariant approach has recently allowed an extension to higher dimensional black 
holes [13, 14], including the tensor modes that only appear in higher dimensions [19]. 

The potential function V for the vector/axial/Regge- Wheeler perturbation is [15, 13, 17] 



v v = f£>\*a + «i + i) 

and for the scalar/polar/Zerilli perturbation the potential is [16, 13, 17] 



V S 



-9Hr + 9m^r - 3m - + 2m z + m" 



(6) 



(7) 



(mr — 3/x) : 

with m = (I — + 2) and I = 2, 3,4, ... is the angular momentum. Note that the I = 0, 1 
modes are special. The former is spherically symmetric and hence, by Birkhoff's theorem, 
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does not describe gravitational radiation. Rather it corresponds to a perturbation changing 
the mass parameter M. The I = 1 mode is in fact pure gauge and arises from a translation. 
We discuss this further in the conclusion below. 

2.3 Boundary conditions and self-adjoint extensions 

Because of the naked singularity the spacetime fails to be globally hyperbolic. In general it 
is far from obvious how to define the dynamics of any field in such a non-globally hyperbolic 
spacetime. However, since the background spacetime we consider is static we can in fact 
define sensible dynamics for linear perturbations. A general prescription for doing so is as 
follows [20, 9, 10, 21]. 

Let A be the spatial derivative part of (5) in a non-globally hyperbolic, static spacetime 

d 2 

A = --^ + V, 0<r*<oo. (8) 

One may view A as an operator acting on the Hilbert space TL = L (r*, dr*) on a timeslice 
Tit orthogonal to the static Killing field. In order to define unitary dynamics, we need a 
self-adjoint extension Ae of A. We will see shortly that choosing a self-adjoint extension 
corresponds in our case to choosing boundary conditions at the naked singularity r* = 0. 
Given a self-adjoint extension, then the time evolution of $ with normalisable initial data 
($0, $o) on S is given by 

&t = cos(4 /2 t)$ + A~ l/2 sin(^ /2 t)$ , (9) 

where $0 = d$/dt\z . One rigorous result that may be proven [20] is that whenever the 
initial conditions are smooth and with compact support, (<I>o, <&o) £ Co°(^o) x Co°(£o), t nen 
&t is smooth everywhere in the spacetime and furthermore within the domain of dependence 
of the initial surface So, <&t agrees with the solution to eq. (5) determined from the initial 
data ($o>^o)- Furthermore, it was shown in [21] that under certain reasonable conditions 
the prescription (9) is the unique way of defining the dynamics. 

Our main interest is in the positivity of Ae or lack thereof. If Ae is positive then the 

1 /2 — 1 /2 1 /2 

dynamics is classically stable since cos(A^ t) and A E sm{A^ t) in (9) become bounded 
operators. Therefore the time evolution of remains bounded for all time and the naked 
singularity is stable at the linearised level. On the other hand, if all possible self-adjoint 
extensions are not positive, then the spacetime is unavoidably unstable. 

Given a self-adjoint extension Ae, one knows that any vector in the Hilbert space may be 
expressed in terms of a basis of eigenvectors of Ae- In particular we could prepare smooth 
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initial data with compact support if we like, even though the individual eigenvectors of Ae 
will satisfy neither of these properties. We therefore turn to a mode analysis of the stability. 
A rigorous proof of stability follows from establishing positivity of Ae [22]. 
More concretely, consider a mode 

$ t = e -^*$(r,) . (10) 

The spatial part then satisfies the Schrodinger equation 

A<$> = -^ + V$> = u 2 <$>. (11) 

Once we fix the boundary conditions at the singularity, the corresponding self-adjoint op- 
erator Ae will be positive if all normalisable solutions to (11) have to 2 > 0, and hence ui 
real. 

Let us determine the possible boundary conditions. Near the singularity, the operator 
for the vector perturbations takes the form 

^v—-^ + ^ + ---, as rwO. (12) 

Therefore, the general solution to (11) behaves as 

$ ~ ai (r; 1/2 + ■■■) + h(r 3 J 2 + •••), as rwO. (13) 

It is immediate to see that both normalisability and self-adjointness require a± = 0. Thus 
there is a unique self-adjoint extension Ae of A in this case, which is defined on the restricted 

1/2 

set of functions that satisfy r* <3>| rit= o = 0. 

The operator for scalar perturbations takes the form 

As ~-^-h + -^ " (14) 
In this case the general solution to (11) near the singularity is 

$ ~ a (rl /2 log- + ■■■)+ b (rl /2 + •••), as n^Q. (15) 
A* 

In this case any choice of (oo, bo) ^ (0, 0) gives normalisable functions and the corresponding 
extension of A is always self-adjoint. From the linearity of the equations involved it is clear 
that (ao,&o) defines the same self-adjoint extension as A(ao,&o)- This equivalence relation 
implies that the family of self-adjoint operators is parameterised by 

q = h G Kpi. (16) 
a 

The positivity of Ae will depend on this parameter. A result of this work will be to show 
for which values of q the singularity is stable. 
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2.4 A critical boundary condition for stability 

We immediately see that the potential Vy for vector perturbations (6) is positive definite. 
Therefore the lowest eigenvalue of Ay in (11) is positive and the spacetime is stable under 
vector perturbations. 

For scalar perturbations the situation is more complicated. Our strategy is as follows. 
Firstly we will show that there is a unique boundary condition, q = qc, such that the 
minimum eigenvalue of As is precisely zero, uj 2 = 0. This critical boundary condition 
separates positive and non-positive self-adjoint extensions. We will then exhibit boundary 
conditions with positive and non-positive spectra on either side of qc- 

The most dangerous mode is the I = 2 mode, as modes with higher angular momentum 
are more positive. Figure 4 shows a plot of the I = 2 potential. For simplicity we will only 
consider the I = 2 case in this subsection. 



400- 




-400- 

Figure 4. The potential Vs(r*) for the I = 2 scalar mode and \i = 1 black hole. The 
potential is unbounded below at the singularity and blows up at an interior point. 

A curious feature of the scalar potential is that the infinite centrifugal barrier appears 
to have been shifted away from the origin to a finite radius rc = 3/i/4. This divergence is 
not a physical effect, but rather an artifact of the variables used to put the perturbation 
equations in Schrodinger form. We will give explicit transformations below to curvature 
variables. The divergence will have the effect of localising the ground state, and hence any 
potential instabilities, near the singularity. This is reasonable given that the perturbations 
are required to be normalisable and that fi is the only length scale in the system. One can 
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check that near the singular point the solutions behave as 

$ ~ Ci[(r - r c f + ■■■]+ C 2 [(r - re)' 1 + •■ ■] as r -► r c , (17) 

where Ci and C2 are constants. The solution that diverges as r — > rc is not normalisable, 
so one must impose that C2 = 0, that is 

<I> ~ (r — rc) 2 as r — > rc* . (18) 

This implies that $ = = at r = rc*- Therefore we may truncate the solution at r = rc 
and match it onto any other solution for r > rc that also satisfies $ = <J>' = at this point. 
In particular, we may take $ = for r > rc- This will be the ground state because V5 is 
positive for r > rc- 

Remarkably, following observations in [14], one may explicitly solve the scalar pertur- 
bation Schrddinger equation in the case when uj 2 = 0. One finds 

r(3fi 3 — 6r 2 |U + 4r 3 ) 



$ = a 



3- 



+ c 4 



/i 3 (4r - 3/i) 

r(13// 3 - 24,u 2 r + 12 / ur 2 ) r(3^ 3 - 6r 2 /i + 4r 3 ) r + //' 



(19) 



3^ 3 (4r-3^) ^i 3 (4r-3//) to r 

where C3 and C4 are constants. We must take the linear combination that does not diverge 
at r = rc = 3/i/4. The well-behaved solution has 



C3 = C4 



7 4 

log 

S 3 9 



(20) 



This solution is seen to have no nodes for r < rc- The expression (19) diverges as r — ► 00. 
However, as indicated above we can take $ = beyond r = rc- Thus we have a normalisable 
uj 2 = ground state. The limit as r — > of this solution therefore gives us the critical 
boundary condition. Comparing with (15) we can read off 

„ , 98 

q c = 2 + logy (21) 

We have at present no conceptual understanding of why there is a marginal mode and why it 
occurs at this value of q. The existence of a boundary condition giving uj 2 = suggests that 
q < qc wm gi ve ground states with co 2 > and hence stable spacetimes and that q > qc will 
give unstable spacetimes. To show this it is sufficient to exhibit one boundary condition with 
q < qc that gives a positive self-adjoint extension and one negative self-adjoint extension 
with q > qc- 

In fact the situation is a little more complicated because q takes values in MP 1 = S 1 
rather than R. There must therefore be at least one other critical boundary condition at 
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which the minimum eigenvalue changes discontinuously. Figure 5 illustrates what happens 
in the simplest case of one discontinous point. This case occurs for perturbations in Anti-de 
Sitter space [12]. 



Unstable 




Figure 5. Stable and unstable boundary conditions as a function of q. 



2.5 A boundary condition for stability 

In this section we show that the boundary condition do = in (15), corresponding to 
q = ±oo, gives a positive self-adjoint extension. The method will be to find a first order 
operator of the form 

« S A + S , (22) 

with S being some smooth function of r. In terms of this S-deformed operator, we formally 
obtain 

($,A$) L 2 = * + / dr* (\D<S>\ 2 + V\<&\ 2 ) , (23) 

L Jr*=0 J V / 

where 

V = V + f^--S 2 , (24) 

where / is the metric function (3) and V is the potential in (8). The expression is formal 
because we should check that the boundary term and integral are not both infinite. 

For any smooth function $ of compact support the boundary term vanishes and inte- 
gration in the second term is finite. If V is now shown to be positive for some appropriately 
chosen S, then the symmetric operator As with domain C£°(r*) is positive definite. Indeed 
we can find a function S that makes manifest the positivity of the symmetric operator A 
for scalar perturbations. We take S to be 

S = --. (25) 
r 



Then we have 



- f(r)U(r) 

Vs = 16(mr - 3 M ) 2 ' {W) 
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where 

fj(r) = 16m (mr + -m 2 + 2m , (27) 

3r z 3 

and where as before m = (I — l)(l + 2). It is clear that Vs is positive definite. 

Note that Vs is unbounded above at r = 3/x/m. For / = 2 we have m = 4 and hence 
this is the same divergence in the potential as we found before. As we saw in the previous 
subsection, states <3? must vanish at this point in order to be square integrable. In this way, 
we see that As is a positive symmetric operator with domain consisting of smooth functions 
of compact support satisfying the regularity condition (18) as r — ► rc- 

Next we need to extend the domain so that As becomes self-adjoint. This is possible 
because a positive symmetric operator always has at least one positive self-adjoint extension, 
known as the Friedrichs extension [20]. We can be completely explicit in our case: Any self- 
adjoint extension corresponding to boundary conditions satisfying 



> , (28) 

r»=0 



will be positive. From (15) we can see that the boundary term is zero if ao = and diverges 
otherwise. Thus we have shown that the boundary condition ao = results in a stable 
spacetime. This boundary condition corresponds to the Friedrichs extension, and may be 
viewed as the simultaneous imposition of generalised Neumann and Dirichlet boundary 
conditions in the sense of [12]. 

2.6 Some boundary conditions for instability 

We may use the Rayleigh-Ritz trial function method to exhibit explicitly a range of values 
for q resulting in unstable spacetimes. Adapted to our context, one has that for any trial 
function x G H we have 

" mi , < -j^y , (29) 

where (,) denotes the inner product of L 2 (r*,dr*). We consider the following set of func- 
tions, depending on three parameters P, Q, S 



X 



3 



(r + Sr 2 ) log - + Pr + Qr 2 
I 1 



for r<^-, (30) 



and vanishing for r > 3/4/i. In this ansatz, the boundary condition at the singularity is 
specified by P. Substituting the ansatz into (29) we minimise the functional with respect 
to Q and S at fixed P. Using these minimum values we see that (%, Asx) is negative for 
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P in the range [1.850, 5.087]. Translating into values of ao using the relation (4) between r 
and it follows that 

q £ [4.393,10.867] Unstable. (31) 

Note that the critical value we found previously has numerical value qc = 4.389, so the trial 
functions give a good approximation to the exact critical value. 



2.7 Energetics 

In this subsection and the following we will relate the variable $ for the scalar perturbations 
to curvature perturbations and to an energy integral. We hope in this manner to clarify 
the presence of an unphysical divergence in the Schrodinger potential and to see if there is 
any sense in which the ao = boundary condition is preferred. 

One may immediately write down an energy which is always finite and conserved and 
furthermore always positive for positive self-adjoint extensions: 



(32) 



where we use (,) to denote the inner product of L 2 (r*, dr*) and the dot denotes a time 
derivative. However, in order to connect this expression to the usual energies that are 
considered in physics one should integrate by parts to obtain a term quadratic in <9 r „<l? 
rather than <J>d r »r« It is at this point that the different boundary conditions produce 
different behaviours. Let us see what happens if we consider energies with terms quadratic 



in d rt ,&. 



There are at least two notions of energy that we could consider. The first is the energy 
associated with the Schrodinger equation 



(33) 



j dr* 




d<S> 


2 




2 




dt 


+ 


<9r* 


+ V S \<$>\ 2 



Use <?i to denote the integrand of the previous expression 

E\ = J £idr* . 



(34) 



By finding a series expansion for the solutions of the Schrodinger equation (11) one finds 
that if ao = the energy is finite whilst if ao / then £\ diverges at the singularity 

|2 



£i 



\ a o\ 



log h 



as 



0. 



(35) 
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The corresponding total energy is therefore infinite. This result might be taken to suggest 
that ao = is the more physical boundary condition. The difference between E$ and E\ 
may be understood as arising upon an integration by parts 

So = S i - ^r$*D$ , (36) 

where D is given in (22) and So is defined analogously to (34). The boundary term is the 
same that arose in (23) and diverges unless ao = 0. 

So far the energies we have considered have been somewhat abstract. However, we can 
clarify the relation of (33) to a more commonly considered definition of energy in general 
relativity. An expression for the energy of a perturbation about a background spacetime is 
given by an integral over a spatial hypersurface 

^ = -^/Gi 2 fc ) n^S (3) , (37) 

(2) 

where G ab is the quadratic variation of the Einstein tensor under the perturbation, re is a 
unit vector orthogonal to the hypersurface, n t = l// 1 / 2 , and ^ is the timelike Killing vector 
£* = 1. The linear variation will vanish because the perturbation satisfies the equation of 
motion. 

For the energy to be conserved, there must be no flux of energy into or out of the 
singularity. Assuming that the perturbation dies off sufficiently rapidly at infinity, the 
change in the energy with time is seen to be 

E 2 = lim i- r 2 f{r) l l 2 G$do {2) . ( 38 ) 
Using, for example, the formalism developed by Chandrasekhar in [17, 18], one can 

(2) 

calculate G b in terms of the scalar metric perturbations. The linearised equations of 
motion for the perturbations are equivalent to the scalar Schrddinger equation. We will not 
review Chandrasekhar's formalism here, although we have used it extensively to obtain the 
results in this subsection and the next. Considering the r — > limit one finds two possible 
behaviours corresponding to the two possibilities in (15). The boundary term (38) may be 
shown to vanish in both cases, so that 

E 2 = 0. (39) 

Thus (37) provides a conserved energy for the perturbations. If we evaluate (37) we again 
find that the different boundary conditions at the singularity give qualitatively different 
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results. The ao = boundary condition has finite energy, whilst the other cases may be 
shown to give 

<?2~ Kl 2 — + ••• as r*^0. (40) 

r* 

Here £2 is defined analogously to (34). Again, the total energy diverges unless ao = 0. Thus 
we see that finiteness of the energy is independent of whether we use E\ or Ei- However, the 
degree of divergence itself is slightly different. We might therefore ask whether there is any 
relation between these two expressions for the energy. Indeed there is. Using, for instance, 
Chandrasekhar's formalism [17, 18] the two expressions may be shown to be identical up to 
total derivative terms 

d 

£2 oc Si + — X 

=> Ei oc E 2 + [X}™ . (41) 

If ao 7^ then 

X~ |a | 2 log 2 — + ••• as r*^0. (42) 
A* 

Therefore the two energies are related by a divergent boundary term. 

2.8 Curvature 

2.8.1 Curvature scalars 

Given that the background is Ricci flat, the natural curvature scalar to consider is the Weyl 
tensor squared. The background has 



C abcd C abcd = ^f. (43) 



It is possible to derive an expression for the linearised perturbation to this curvature scalar 
in terms of the Schrodinger variable $ using, for example [17, 18]. One finds 



S (c abcd C abcd ) (r) = 



— 12/xm (m + 2) 



$(r) (6r/x + V)/(r)V 2 f <S>(r')dr' 



r 6 (m + 2)r 8 J f{r') l l 2 (mr'-^) 



S 



e- iujt Pi(cos9), 
(44) 



where we have restored the angular and time dependence and without loss of generality we 
have restricted to axisymmetric perturbations for simplicity. 

There are at least two observations to make about equation (44) . Firstly, we can see that 
there is no physical divergence at the point r = 3/x/m, at which the Schrodinger potential 
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diverges. Given that $ must go to zero at this point, the integral in (44) is finite. Secondly, 
taking the limit r->0we find 



6 (c ahcd C abcd ) ~ a fllog- + . 



+ b o ( 4 + • • • 



as 



0. 



(45) 

/i / \r u / 
Therefore we see that if ao = the perturbed curvature has the same degree of divergence 
as the background, whilst if ao 7^ then the divergence of the perturbed curvature is greater 
than that of the background. 



2.8.2 Weyl scalar s 

Perturbations about black hole spacetime are often considered using the formalism of Weyl 
scalars. This is particularly useful for the case of rotating black holes where the Newman- 
Penrose formalism allows a separation of variables of the perturbation equations [24, 25]. 
In this subsection we briefly recast our analysis in terms of a Weyl scalar. 

An explicit transformation is known [17, 18] between the variable $ satisfying the scalar 
Schrodinger equation (the Zerilli equation) and the perturbed Weyl scalar 
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In the Schwarzschild background the only nonvanishing Weyl scalar is 
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The relationship between the perturbed Weyl scalar and the scalar perturbation, including 
the angular and time dependence, is [17, 18] 
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This equation was derived in [17, 18] by working in a specific gauge, but given that 5^q 
is gauge invariant for the Schwarzschild background, the relation should also hold for our 
gauge invariant <E>. 
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In (49) we can see the appearance of a divergence at the finite radius r = 3/x/m. For 
the Weyl scalar to remain bounded at this radius, <3? must vanish at r = 3fi/m. This is 
precisely the condition that we found previously due to a divergence in the scalar potential. 
Therefore, we see explicitly that the unphysical divergence in the potential arises in the 
change of variables relating the metric perturbation variable to physical curvatures. 

If we consider the r — > limit of $ (15) then we find that the Weyl scalar perturbation 
behaves as 



3 Discussion 

The linearised dynamics of gravitational perturbations about the negative mass Schwarzschild 
spacetime are given a well defined dynamics by specifying boundary conditions at the sin- 
gularity. There is a one parameter choice of possible boundary conditions corresponding 
to self-adjoint extensions of the Hamiltonian. We considered the stability of the spacetime 
as a function of the boundary condition. We have shown that there is a critical boundary 
condition separating stable and unstable spacetimes. 

Amongst the possible boundary conditions, we have seen that there is one particular 
choice, ao = in (15), which gives perturbations with finite physical energies E\ and 
E2. Furthermore, these perturbations induce curvature perturbations with the same degree 
of divergence at the singularity as the background curvature. All the other boundary 
conditions result in perturbations with infinite energy which have an enhanced curvature 
divergence at the singularity. Therefore, the ao = boundary seems to be more physical 
than the others. We saw that this boundary condition gave a stable spacetime. 

The main conclusion following from these results is that the negative mass Schwarzschild 
phantom can be perturbatively stable and in particular is stable with what appears to be 
the most physical choice of boundary conditions. This is perhaps counter to intuition, given 
the existence of nonlinear instabilities that we reviewed in the introduction. 

There are two immediate limitations of our calculations. Firstly, very near to the sin- 
gularity it is likely that a linearised analysis of perturbations is not valid. This is because 
higher order terms in the equations of motion will contain curvature tensors that diverge as 
r* — ► 0. By making the perturbations sufficiently small one might hope that the linearised 
approximation will be valid down to very small radii. Perhaps in this case the nonlinear 
effects near the singularity could be absorbed into an effective boundary condition. Indeed 




(50) 
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it would be interesting if a preferred boundary condition is selected in this way. 

Secondly, near the singularity the curvatures become large and it seems likely that 
quantum corrections to the spacetime will be important. Again, such corrections could 
perhaps be renormalised into an effective boundary condition. In this connection, it is 
interesting to note that there are arguments suggesting that the negative mass Schwarzschild 
singularity should not be resolved by quantum corrections [23]. 

We would now like to return to the issue of the I = 1 scalar mode. As mentioned above, 
this is a 'translational zero mode' which could be interpreted as allowing the singularity to 
move from its initial position. Because the mass of the singularity is negative, one expects 
such a motion to carry negative kinetic energy. This could clearly lead to various nonlinear 
instabilities including that envisioned by Bondi, which we reviewed in the introduction. 

There are various possible extensions of our calculations. Adapting to higher dimensional 
spacetimes should be straightforward, the necessary formalism may be found in [13, 14, 
19]. Back in four dimensions it may be possible to extend the calculations to the case of 
rotating Kerr black holes. Remarkably, the linearised stability of positive mass Kerr black 
holes is a tractable problem and stability has been proven [24, 25, 18, 26]. Negative mass 
rotating black holes would presumably also have a negative moment of inertia and therefore 
potentially suffer from additional 'spin up' instabilities. 

More speculatively, there is an intriguing similarity between the behaviour of the per- 
turbations we have considered near the singularity and the asymptotic behaviour of pertur- 
bations about Anti-de Sitter space [12]. Perhaps there is a version of 'holographic renormal- 
isation' [27, 28, 29] which applies to the naked singularity and allows one to absorb some 
divergences arising from integrations near the singularity? 
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